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I.  INTRODUCTION 


In  the  post  Gulf  War  era,  significant  effort  has  been  expended  in  the  area  of  Non-Coopera¬ 
tive  Target  Identification  (NCTI).  The  techniques  proposed  for  use  with  radar  systems  all  have  a 
common  feature:  the  need  to  use  target  aspect  angle  with  respect  to  (w.r.t.)  the  direction  of  arrival 
of  the  radar  signal  to  reduce  the  problem  to  a  manageable  size.  This  report  will  present  a  rigorous 
development  of  the  mathematics  for  solving  this  problem  exactly  given  target  position  and  atti¬ 
tude  data;  and  then,  develop  the  mathematics  for  solving  this  problem  when  only  the  radar’s  posi¬ 
tional  measurement  of  the  target  is  available.  Finally,  data  will  be  presented  to  show  the  practical 
accuracy  limits  of  aspect  angle  calculations  using  radar  data. 

n.  DEFINITIONS 

The  basic  problem  is  to  find,  w.r.t.  the  aircraft,  the  two  angles  that  describe  the  direction  of 
arrival  of  a  radar  beam  at  the  aircraft.  Call  these  angles  0  and  (|)  where  0  is  the  angle  between  the 
aircraft’s  nose  and  the  projection  of  the  radar  beam  into  the  aircraft’s  wing  plane  measured  in  the 
counter-clockwise  direction  when  viewed  fi'om  above  and  (j)  is  the  angle  between  the  aircraft’s 
wing  plane  and  the  vector  that  describes  the  radar  beam  striking  the  aircraft  measured  in  the 
clockwise  direction  when  viewed  from  the  port  wing. 

The  information  given  to  solve  this  problem  is  the  radar’s  measurement  of  the  aircraft’s  cur¬ 
rent  position  in  a  Rng,  Az,  El  coordinate  system  referenced  to  the  radar.  Call  this  point  P  and 
denote  it  as  a  function  of  Rng,  Az,  and  El  by  P  (Rng,  Az,  El).  Also  given  is  the  aircraft’s  spacial 
orientation  w.r.t.  the  radar  in  the  form  of  a  heading  angle,  a  pitch  angle,  ti;  and  a  roll  or  bank 

angle,  The  heading  angle  is  the  angle  between  true  North  and  the  aircraft’s  centerline  in  the 
horizontal  plane  of  the  radar  measured  in  the  clockwise  direction.  The  pitch  angle  is  the  angle 
between  the  radar’s  horizontal  plane  and  the  aircraft’s  center  line  measured  in  the  clockwise 
direction.  The  roll  angle  is  the  angle  between  the  radar’s  horizontal  plane  and  the  aircraft’s  port 
wing  centerline  measured  in  the  clockwise  direction. 
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ni.  CLOSED  FORM  SOLUTION 


A.  6  and  (|>  w.r.t.  the  Aircraft’s  Reference  Frame 

As  a  first  step  in  solving  this  problem,  assume  a  right-handed  Cartesian  coordinate  sys¬ 
tem  with  origin  at  the  aircraft  and  unit  vectors  ila ,  ,  and  Za  where  ^a  points  in  the  direction  of 

the  aircraft’s  nose,  ya  points  out  the  aircraft’s  port  wing,  and  Za  points  up  from  the  y^  plane  in 

f 

the  right-handed  sense  (Fig.  1).  Let  'V'^ar  be  the  pointing  vector  from  the  aircraft  to  the  radar  and 

be  expressed  in  the  form  Var  =  Varx  Xa  +  V^Ry  ya  +  Varz  Za  •  Then  from  Figure  1,  the  aspect 
angles  6  and  <(>  are: 


0  =  cos 


-1 


ARx 


)'  +  (Var/. 


(1) 


({)  =:  tan 


-1 


^ARz 


[J(^. 


f+(VARyf. 


(2) 


Figure  1.  Aircraft  Aspect  Angles 
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Thus,  the  major  step  in  the  determination  of  the  direction  of  arrival  of  the  radar  beam 
w.r.t.  the  aircraft  is  the  transformation  of  the  radar’s  measurement  of  the  aircraft’s  position  into  a 
pointing  vector  from  the  aircraft  to  the  radar  in  the  coordinate  space  described  above. 

B.  Development  of  the  Radar  to  Aircraft  Pointing  Vector,  Vra  ,  w.r.t  the  Radar’s 

Reference  Frame 

The  position  data  from  the  radar  is  typically  in  a  Rng,  Az,  El  coordinate  system  refer¬ 
enced  to  the  radar.  As  the  Rng,  Az,  El  coordinate  system  is  a  left-handed  system,  three  transfor¬ 
mations  of  the  radar  data  point  P  (Rng,  Az,  El)  will  be  required  to  map  it  into  the  aircraft’s 
coordinate  system. 

The  first  transformation  is  to  map  the  Rng,  Az,  El  data  into  a  right-handed  East,  North, 
Up  coordinate  system  referenced  to  the  radar.  This  is  a  standard  transformation  in  the  radar  field 
and  is  given  by: 

East  or  X  =  Rng  sin  (Az)  cos  (El)  (3a) 

North  or  y  =  Rng  cos  (Az)  cos  (El)  (3b) 

Up  or  z  =  Rng  sin  (El) .  (3c) 


From  this,  a  pointing  vector  from  the  radar  to  the  aircraft,  Vra,  can  be  defined  as: 

Vra  =  sin(Az)cos(El)x  +  cos(Az)cos(El)y  +  sin(El)z.  (4) 

C.  Development  of  the  Aircraft  to  Radar  Pointing  Vector,  Var  ,  w.r.t.  the  Aircraft’s 
Frame  of  Reference 

The  next  transformation  is  to  transform  the  vector,  Vra,  into  a  pointing  vector,  Var  , 
from  the  aircraft  to  the  radar  in  an  East,  North,  Up  coordinate  system  centered  at  the  aircraft. 

Since  the  aircraft’s  position  is  also  the  head  of,  Vra,  Var  will  simply  be  the  additive  inverse  of 
Vra  or 


Var  =  -1  *  Vra. 


(5) 


The  final  transformation  is  to  transform  the  vector  Var  into  the  vector  Var  as  defined 
earlier.  While  this  is  not  an  extremely  difficult  transformation,  it  is  rather  involved  since  it  is 
effectively  a  three  axis  of  rotation  rigid  body  problem. 
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Using  vector  analysis,  this  type  of  transformation  can  be  written  as: 


(6) 


where: 

X,  y,  z  are  the  axes  labels  for  the  original  coordinate  system, 
x',  y',  z'  are  the  axes  labels  for  the  rotated  coordinate  system, 

and  J  is  the  transformation  matrix  whose  columns  consist  of  the  direction  cosines  of  the  unit 
vectors  x' ,  y' ,  z'  w.r.t.  the  old  coordinate  system. 

Now  the  only  remaining  piece  of  the  puzzle  is  the  derivation  of  the  transformation  matrix  J. 

The  transformation  matrix  J  can  be  derived  through  three  successive  transforms  since 
each  rotation  is  a  linear  operation  and  therefore  separable.  The  first  transform  is  the  result  of  the 
heading  angle,  |,  and  is  a  clockwise  rotation  of  the  original  x,  y,  z  coordinate  system  about  the 
z-axis.  This  is  shown  in  Figure  2.  An  adjustment  must  be  made  to  |  though  due  to  a  reference  dif¬ 
ference  between  the  radar  and  aircraft  coordinate  systems.  The  angle  |  is  referenced  to  North,  but 
the  x-axis  of  the  radar’s  coordinate  system  is  aligned  to  East.  This  results  in  the  need  for  a  new 
heading  angle,  where  ^  -  90  degrees. 


Figure  2.  Rotation  due  to  ^ 
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Using  Figure  2,  the  transformation  is: 


X 

y 

z 


cos(l') 

-sin(^’) 


0 


sin(^')  0 
cos (5')  0 
0  1 


(7) 


The  next  transform  is  the  result  of  the  pitch  angle,  t],  and  is  a  clockwise  rotation  of  the 
x',  y',  z'  coordinate  system  about  the  y'-axis.  This  is  shown  in  Figure  3. 


Figure  3.  Rotation  due  to  r| 
Using  Figure  3,  the  transformation  is: 


X’ 

cos(t])  0 

-sin(r|) 

x" 

y' 

= 

0  1 

0 

y" 

X 

_sin(r|)  0 

cos(t]) 

_z"_ 
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The  final  transform  is  the  result  of  the  roll  angle,  ^  and  is  a  clockwise  rotation  of  the  x", 
y",  z"  coordinate  system  about  the  x"  axis.  This  is  shown  in  Figure  4. 


Using  Figure  4,  the  transformation  is: 


x" 

0  0  0 

y" 

= 

0  cos(^)  sin(^) 

z" 

_1  -sin(C)  cos(^) 

(9) 


Now  substituting  Equations  (8)  and  (9)  into  (7)  yields: 

cos(§')cos(t])  cos(^')sin(r|)sin(^)+sin(|')cos(^) 

-sin  ( I' )  cos  (ti  )  -sin  (!' )  sin  (“n )  sin  ( t) + cos  ( )  cos  (^) 
sin(Ti)  -cos(Ti)sin(C) 

-cos  (!' )  sin  (q )  cos  (^) + sin  (!')  sin  (^) 

sin  (!')  sin  (il )  cos  (C) + cos  (^')  sin(^) 
cos(t])cos(^) 

Equation  (10)  provides  the  transformation  from  the  aircraft’s  coordinate  system  to  the 
radar’s  coordinate  system  which  is  the  inverse  of  the  desired  transformation.  Multiplying  both 
sides  by  J  ‘  yields: 


"'a 

Ya 

Za 


(10) 
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(11) 


Normally,  computing  the  inverse  of  a  matrix  with  complicated  terms  like  the  ones  in  J 
would  be  quite  involved;  however,  since  this  is  a  linear  orthogonal  transform,  the  inverse  of  the 
transformation  matrix  is  equal  to  the  transpose  of  the  transformation  matrix,  And  corre¬ 

spondingly: 


cos(^')cos(ti) 

cos  (!' )  sin  (ti  )  sin  (^)+ sin  (!')  cos(^) 
-cos(^')sin(ii)cos(^)+sin(|')sin(^) 


-sin(|')cos(Ti)  sin(q) 

-sin(i')sin(Ti)sin(^)+cos(|')cos(^)  -cos(Ti)sin(^) 
sin(|')sin(ii)cos(^)+cos(^')sin(^)  cos('q)cos(^) 


(12) 


This  completes  all  the  pieces  needed  to  transform  Vra  into  Var  .  This  transform  can 
be  written  as: 


Var  =  J'^-Vra)  (13) 

Which,  when  substituting  in  Equation  (4)  and  expanding,  yields: 

Varx  =  -cos(?')cos(Ti)sin(Az)cos(El)  + sin(^')cos(q)cos(Az)cos(El) 

-sin(Ti)sin(El)  (14a) 

^ARy  =  -(cos(^')sin(ti)sin(^)  +  sin(|')cos(^))sin(Az)cos(El) 
-(-sin(§')sin(Ti)sin(^)  +  cos(^')cos(^))cos(Az)cos(El) 

+  cos(T|)sin(^)sin(El)  (14b) 
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Varz  =  -(-cos(l')sin(Ti)cos(^) +  sm(i')sm(^))sin(Az)cos(El) 
-(sin(|')sm('n)cos(U  +  cos(|')sin(^))cos(Az)cos(El) 
-cos(Ti)cos(^)sm(El) .  (14c) 

These  equations  can  then  be  substituted  into  Equations  (1)  and  (2)  for  the  final  aspect 
angles  calculation. 

rV.  PRACTICAL  CONSIDERATIONS  FOR  ASPECT  ANGLES  CALCULATION 

In  the  above  derivation,  complete  knowledge  of  the  orientation  vector  of  the  aircraft  was 
assumed.  For  a  well  instrumented  test  aircraft,  this  is  a  valid  assumption;  unfortunately,  it  breaks 
down  in  a  combat  situation.  In  this  case,  the  aircraft’s  orientation  vector  must  be  calculated  from 
the  radar’s  positional  measurements. 

The  need  to  calculate  an  aircraft’s  complete  orientation  vector  from  radar  positional  mea¬ 
surements  presents  some  difficult  problems  in  the  roll  component.  Calculating  the  aircraft  s  roll 
angle  would  require  a  quadratic  curve  fitter,  good  estimates  of  the  aircraft’s  velocity,  and  turn 
radius  versus  roll  angle  curves  for  all  classes  of  aircraft.  Coupling  this  computational  complexity 
with  the  fact  that  most  high  roll  maneuvers  occur  in  an  end  game  scenario  when  weapons  have 
already  been  committed,  leads  to  the  conclusion  that  in  an  application  sense  it  is  useful  and  rea¬ 
sonable  to  set  the  roll  component  to  zero. 

Once  this  simplification  has  been  made  the  heading  and  pitch  components  of  the  orientation 
vector  can  be  calculated  from  simple  linear  travel.  Let  Pc(Rng,  Az,  El)  be  the  radar  s  current  mea¬ 
surement  of  the  aircraft’s  position  w.r.t.  radar,  and  let  Pp(Rng,  Az,  El)  be  the  radar’s  previous 
measurement  of  the  aircraft’s  position.  Then,  using  Equation  (3),  the  aircraft’s  motion  can  be 


described  by: 

Ax  =  Xj.-Xp 

Ay  =  Yc-Yp 
Az  =  Zp-Zp. 

From  this,  the  pitch  and  heading  can  be  calculated  as: 
^  =  tan~^(Ay/Ax) 


^  (  Az/V(Ax)^  +  (Ay)^  I 
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It  is  worth  noting  that  error  components,  unrelated  to  roll,  do  exist  in  these  calculations.  The 
velocity  based  heading  angle  will  differ  from  the  actual  heading  angle  by  an  amount  known  as  the 
crab  angle.  The  crab  angle  is  defined  as  the  difference  between  an  aircraft’s  direction  of  travel  and 
the  projection  of  its  centerline  and  is  typically  caused  by  wind  velocity.  This  angle  is  small  at  high 
aircraft  velocities  but  can  induce  errors  at  lower  velocities.  In  light  of  the  large  variation  in  wind 
velocities  at  various  altitudes,  this  error  cannot  be  readily  predicted.  Also,  the  actual  pitch  angle 
will  vary  from  the  calculated  climb  angle  by  an  amount  known  as  the  Angle-of-Attack  (AOA). 
The  AOA  results  from  the  need  for  Jet  aircraft  to  pitch  up  in  order  to  maintain  flight.  This  angle 
varies  with  velocity  and  airframe  but  could  be  predicted  once  the  data  is  provided  to  an  NCl'l 
algorithm  and  the  airframe  is  estimated  in  the  algorithm.  In  general  though,  the  use  of  climb  angle 
for  pitch  angle  is  sufficient. 

The  primary  advantage  of  setting  the  roll  component  to  zero  is  this  allows  a  major  simplifi¬ 
cation  of  the  aspect  angle  problem.  Under  this  condition,  the  problem  becomes  one  of  simple 
geometry  and  trigonometry  as  can  be  seen  in  Figures  5  and  6.  Using  these  figures,  the  new  equa¬ 
tions  for  6  and  (|>  are  seen  to  be: 

e  -  p  - 1  (18) 

(|) »  -1(ti  cos(0)  +  El)  (19) 


where: 


p  =  -1  *  Az,  which  is  the  back  azimuth  from  the  aircraft  to  the  radar, 
El  is  the  elevation  angle  of  the  aircraft  w.r.t.  the  radar. 
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Figures.  Simplified  Estimate  of  Q 


Figure  6.  Simplified  Estimate  of^ 

As  will  be  shown  in  the  next  section,  this  method  has  been  found  to  be  within  ±5  degrees  of 
the  actual  value  for  either  aspect  angle.  Given  that  most  NCTI  algorithms  search  over  at  least  this 
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much  angle  space,  this  is  not  an  unreasonable  result.  This  is  especially  true  since  it  is  not  clear  at 
this  time  if  the  roll  component  can  be  calculated  with  sufficient  accuracy  to  improve  these 
bounds. 

V.  AIRCRAFT  TESTS  INVOLVING  ASPECT  ANGLES 

NCn  research  and  development  tests  were  performed  using  the  Multi-Role  Survivable 
Radar  (MRSR)  at  Redstone  Arsenal,  AL.  The  aircraft  involved  were  the  F-15A  and  F-14A. 
MRSR  is  a  highly  accurate,  3-Dimensional  (3-D)  (provides  Rng,  Az,  and  El)  radar  which  pro¬ 
vides  tracking  and  high  range  resolution  data  at  a  2  second  rate.  As  an  input  to  NCTI  algorithms, 
MRSR  must  provide  an  estimate  of  the  target’s  aspect  angles.  In  order  to  evaluate  the  accuracy  of 
this  estimate,  positional  and  attitude  truth  data  had  to  be  collected. 

The  best  available  method  for  obtaining  positional  and  attitude  truth  data  is  the  Advanced 
Range  Data  System  (ARDS).  This  is  a  precise  (encrypted  P-code)  Global  Positioning  System 
(GPS)  receiver  integrated  with  an  Inertial  Reference  Unit  (IRU).  The  system  collects  positional 
and  attitude  data  at  a  10  Hz  rate.  When  used  with  a  reference  receiver  to  provide  differential  cor¬ 
rections,  this  system  can  provide  approximately  1  meter  positional  accuracy  and  approximately 
1  degree  angular  accuracy  in  the  roll,  pitch,  and  yaw  attitude  angles.  When  operated  in  non-differ¬ 
ential  mode,  the  positional  accuracy  degrades  to  approximately  16  meters  with  no  degradation  in 
attitude  accuracy. 

Both  aircraft  in  the  test  were  outfitted  with  an  ARDS  unit  contained  in  an  AIM-9  pod  and 
operated  in  the  non-differential  mode.  MRSR  track  and  ARDS  data  were  collected  on  these  air¬ 
craft.  The  aircraft  were  flown  in  racetrack  patterns  around  the  radar  so  that  data  over  a  wide  range 
of  aspect  angles  were  collected.  Given  the  accuracies  of  the  ARDS  unit,  the  aspect  angles  calcula¬ 
tions  should  be  accurate  to  approximately  1  degree  in  0  and  (j>.  The  error  induced  by  the  positional 
accuracy  of  non-differential,  P-code  GPS  was  insignificant  for  the  target  ranges. 

The  table  shows  the  differences  between  the  aspect  angles  as  calculated  from  MRSR  mea¬ 
sured  data  using  the  equations  in  Section  IV  and  as  calculated  from  ARDS  data  using  the  equa¬ 
tions  in  Section  IE.  The  aspect  angles  differences  were  also  calculated  with  GPS  positional  data 
substituted  for  MRSR  data  in  the  equations  in  Section  IV.  This  gives  the  best  estimate  of  the 
errors  resulting  from  the  assumptions  stated  in  Section  IV  for  determining  aspect  angles  with  only 
positional  information. 

The  table  is  also  broken  down  into  inbound  and  outbound  data.  This  illustrates  the  effect  of 
the  aircraft’s  AOA  above  its  line  of  flight.  Note  that  the  (j)  error  is  positive  (the  underside  of  the 
aircraft  is  visible)  on  inbound  runs  and  negative  (the  topside  is  visible)  on  outbound  runs.  This 
angle  error  is  dependent  on  the  type  of  aircraft,  its  velocity,  and,  obviously,  on  whether  the  air¬ 
craft  is  inbound  or  outbound.  The  velocity  for  these  flights  was  kept  at  a  constant  270  knots  indi¬ 
cated  air  speed.  An  area  for  future  research  will  be  to  attempt  to  develop  a  method  to  reduce  this 
error  given  only  the  aircraft’s  velocity  and  length. 
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Table.  Aspect  Angle  Accuracies 


F-15  aspect  angles  measured  by  MRSR  compared  to  IRU/GPS  aspect  angles 

4)1,  avg  err  =  3.28* 

0™.  avg  err  =  0.52° 

<t)i„  stdev  =  5.14® 

Stdev  =  3.86® 

F-15  aspect  angles  measured  by  P-code  GPS  positional  data  compared  to  IRU/GPS  aspect  angles 

®it,  avg  err  =  -0.45® 

4)i„  avg  err  =  3.26° 

4>m.  avg  err  =  -2.89° 

<});„  stdev  =  2.18® 

4)„„  stdev  =  1.68° 

F-14  aspect  angles  measured  by  MRSR  compared  to 

IRU/GPS  aspect  angles 

4>i,  avg  err  =1.97° 

0™.  avg  err  =  -0.56° 

4>„«  avg  err  =  -4.82° 

(j)j„  Stdev  =  5.46® 

4)„„  stdev  =  5.79° 

F-14  aspect  angles  measured  by  P-code  GPS  positional  data  compared  to  IRU/GPS  aspect  angles 

4)i.avg  err  =  1.71° 

0™.  avg  err  =  -0.27° 

4>„„  avg  err  =  -3.87° 

(j),^  Stdev  =  1.50® 

4>„„  stdev  =  5.57° 

I 
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